DEVELOPING  A  MATHEMATICAL  MODEL  FOR  PREDICTING  TOTAL 
OPERATORS'  WAGES  IN  PUBLIC  TRANSIT 


by 


RAJ  KUMAR  GUPTA 
B,  S.  (M.E.XHons.)  University  of  Anand,  India,  1961 


A  MASTER'S  THESIS 


submitted  in  partial  fulfillment  of  the 


requirements  for  the  degree 


MASTER  OF  SCIENCE 


Department  of  Industrial  Engineering 


Kansas  State  University 
Manhattan,  Kansas 


1965 


Approved 


Major/ Professor 


LD 

2G0f 

T'/ 


^  TABLE  OF  CONTENTS 


SYMBOL  AND  DEFINITIONS   1 

HYPOTHESIS  2 

INTRODUCTION  . 3 

DEVELOPMENT  OF  STRAIGHT  RUNS   8 

CONSTRUCTING  OF  SPLIT  RUNS 13 

SPLIT  RUN  SELECTION 18 

SELECTION  OF  PIECES  OF  WORK  NOT  USED  FOR  SPLIT  RUNS  .  .  20 
ALLOCATING  PENALTIES  TO  THE  UNUSED  PIECES  OF  WORK  ...  21 

DESIGN  OF  EXPERIMENT  .....  22 

MULTIVARIATE  AND  BIVARIATE  CORRELATION  ANALYSIS  ....  25 

TWO-WAY  ANALYSIS  OF  VARIANCE   31 

BIBLIOGRAPHY .38 

ACKNOIVLEDGMENTS 39 

APPENDIX 40 


TABLE  OF  CONTENTS 

SYMBOL  AND  DEFINITIONS   1 

HYPOTHESIS  2 

INTRODUCTION • 3 

DEVELOPMENT  OF  STRAIGHT  RUNS 6 

CONSTRUCTING  OF  SPLIT  RUNS 13 

SPLIT  RUN  SELECTION   18 

SELECTION  OF  PIECES  OF  WORK  NOT  USED  FOR  SPLIT  RUNS  .  .  20 
ALLOCATING  PENALTIES  TO  THE  UNUSED  PIECES  OF  WORK  ...  21 

DESIGN  OF  EXPERIMENT  22 

MULTIVARIATE  AND  BIVARIATE  CORRELATION  ANALYSIS  ....  25 

TIVO-WAY  ANALYSIS  OF  VARIANCE 31 

BIBLIOGRAPHY 38 

ACra^OlVLEDGMENTS 39 

APPENDIX 40 


LIST  OF  FIGURES 

FIGURE     TITLE  PAGE 

1  DAILY  VOLUME  OF  TRAFFIC  CURVE  5 

2  COMPUTER  OUTPUT  -  STRAIGHT  RUN  MAKING  12 

3  COMPUTER  OUTPUT  -  COMBINING  PIECES  OF  WORK  ...  17 

4  ACTUAL  TOTAL  PAY  TIME  FOR  OPERATORS 24 

5  PREDICTED  TOTAL  PAY  TIME  FROM  EQUATION  1   ...  27 

6  PREDICTED  TOTAL  PAY  TIME  FOR  EQUATION  2  .  .  .  .  28 

7  ANALYSIS  OF  VARIANCE  TABLE  34 

8  ACTUAL  TOTAL  PAY  TIME 35 

9  COMPUTATION  FOR  ANALYSIS  OF  VARIANCE  36 


DEFINITIONS 


Straight  Run 
Piece  of  Work 

Split  Run 
Pull  Out  Time 

Pull  In  Time 

Travel  Time 

Platform  Time 
Turn  In 

Guarantee  Time 

Allowance 

Pay  Time 
Spread  Time 

Relief  Point 


A  regular  run  having  no  unpaid  breaks. 

Any  work  assigned  to  an  operator  which  pays  less 

than  eight  hours. 

A  regular  run  having  unpaid  breaks. 

The  starting  time  of  a  piece  of  work,  straight  run, 

or  a  split  run. 

The  finishing  time  of  a  piece  of  work,  straight  run, 

or  a  split  run. 

The  actual  running  time  from  the  garage  to  point  of 

relief. 

The  time  during  which  an  operator  is  in  charge  of  a  bus. 

A  specific  amount  of  time  which  the  operator  receives  for 

every  straight  run  and  each  piece  of  work  in  his  assignment. 

The  time  paid  for  but  not  worked  to  make  up  a  minimum 

guarantee  for  certain  classes  of  work. 

A  time  equal  to  either  the  travel  time,  the  guarantee 

time,  or  both. 

The  total  pay  time  for  a  specific  run. 

The  total  working  time  plus  the  unpaid  time  of  a  split 

run. 

Relief  point  is  a  certain  location  or  locations  on  each 

route  where  operators  may  be  relieved  or  assigned  to  the 

bus. 


FlYPOTliFSIS 

The  task  of  public  transit  companies  is  becoming  immeasurably 
greater  and  more  complex  due  to  the  declining  trend  in  the  usage  of  their 
services.  One  of  the  urgent  needs  of  todny  is  that  of  increasing  the 
efficiency  of  the  transit  services  and  to  make  them  more  attractive  to 
the  public.   In  view  of  the  constant  increase  in  the  cost  of  labor,  there 
is  an  urgent  need  for  reducing  operating  costs.  Proper  scheduling  techni- 
ques yielding  efficient  and  effective  schedules  would  play  a  major  role 
in  reducing  operating  costs.  Hence,  the  aim  of  transit  companies  has 
been  rescheduling  services  as  a  means  of  maintaining  operations  on  a 
profitable  basis  in  face  of  strong  competition  from  automobile  users. 

The  objective  of  this  thesis  is  to  develop  a  mathematical  model  for 
predicting  the  total  pay  time  to  the  drivers  of  public  transit  industry 
in  a  general  functional  forms  and  parameter  values.  It  has  been  found 
that  there  are  various  constraints  and  restrictions  in  the      ir.nnagement 
agreement  on  which  always  a  controversy  is  going  on  between  management  and 
trade  union.  Further,  it  is  felt  necessary  to  explore  the  sensitivity  of 
the  model  to  some  general  functional  equation  in  terms  of  parameter  values. 
The  carrying  out  of  computation  of  the  model  is  done  in  1620  IBM  Computer. 
The  total  pay  time  functions  are  plotted  for  each  combination  of  parameter 
values.  Also  statistically  by  conducting  a  two-way  analysis  of  variance 
test  it  is  intended  to  show  which  of  these  variables  affect  more  signifi- 
cantly the  total  pay  time  function. 
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INTRODUCTION 

The  automobile  vehicle  brought  a  new  concept  to  transportation, 
flexibility.  The  first  twenty  years  of  the  20th  Century  were  devoted 
to  perfecting  the  mechanical  design  of  the  vehicle.  Few  people  recognized 
the  potential  of  the  motor  vehicle  or  anticipated  the  revolution  it  would 
cause  in  American  economy  or  in  cultural  and  social  values.  It  is  un- 
necessary to  enumerate  the  tremendous  changes  this  brought  in  America 
and  the  great  advance  in  economic  progress  it  made  possible.  No  longer 
are  there  isolated  areas  in  the  nation  inaccessible  to  the  amenities  of 
life. 

In  many  cities  the  bus  system  is  one  of  the  major  means  of  trans- 
portation into  or  within  the  city.  There  are  various  reasons  for  adoption 
of  buses  in  transit  operations  among  which  a  few  are  as  follows: 

(1)  Flexibility  in  traffic, 

(2)  Individual  power  supply, 

(3)  Ability  to  pass  each  other, 

(4)  Through  service  and  off  route  locations, 

(5)  Ability  to  combine  routes  with  one  vehicle  rendering  the  service, 

(6)  Low  initial  cost. 

From  the  last  few  years  the  number  of  intracity  bus  passengers  has 
been  steadily  decreasing.  There  are  various  reasons  for  this  decrease, 
but  the  decline  in  the  usage  of  the  buses  and  the  increase  of  automobiles 
on  the  road  have  not  only  created  problems  for  transit  companies,  but  also 
for  the  city  planners.  For  reversing  this  trend  many  solutions  have  been 
suggested,  but  time  is  needed  to  assess  the  validity  of  these  solutions. 
However,  even  those  transit  companies  which  were  making  tremendous  profits 
are  now  facing  a  problem  to  stay  in  the  business  because  of  the  decrease 


in  the  volume  of  bus  passengers  and  increase  of  automobile  users. 

Since  the  transit  companies  are  service  companies,  a  major  portion 
of  their  revenues  goes  as  wages  to  operators.  Thus,  it  has  been  felt  that 
the  reduction  in  operating  cost,  through  rescheduling  is  the  only  way  for 
transit  companies  to  maintain  their  operation  on  a  profitable  basis.  Yet, 
the  best  manual  methods  fall  short  of  solving  the  problem,  because  com- 
putation cannot  be  made  quickly  enough  to  keep  pace  with  the  changing 
conditions.  As  a  result,  high  speed  computers  have  proved  the  best  means 
to  overcome  these  difficulties. 

The  theoretical  analysis  involved  in  scheduling  has  been  completed 
in  previous  research  and  computer  programs  to  carry  on  the  scheduling 
operations  has  been  developed  which  have  been  modified  here  to  include 
each  company's  constraints  and  variables. 

The  number  of  variables  and  restrictions  is  the  main  difficulty 
encountered  in  effectively  programming  a  daily  work  schedule  for  each 
operator.  Operating  variables,  trip  times,  and  restrictions  of  the  union- 
management  agreement  are  a  few  of  the  factors  which  have  to  be  considered 
in  making  an  effective  daily  schedule. 

One  of  the  biggest  problems  faced  by  the  transit  industry  is  the 
variable  nature  of  the  daily  demand  of  their  services.  Traffic  is  at 
the  peak  in  the  mornings  and  late  afternoons  and  tends  to  fall  off  in  the 
early  mornings,  noon,  and  late  evenings  and  on  week  ends.  To  meet  these  peak 
traffic  requirements,  it  becomes  necessary  to  maintain  a  big  fleet  of  buses 
in  excess  of  that  which  otherwise  would  be  needed.  The  graph  in  Figure  1, 
shows  the  requirement  of  the  number  of  buses  (motor  coaches)  needed  during 


^Samy  E.  G.  Elias,  "A  Digital  Computing  Solution  to  the  Transit  Operation 
Assignment  Problems,"  Doctoral  Thesis,  Oklahoma  State  Univ.,  1960, 
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the  day  in  a  hypothetical  case.  From  the  Figure,  it  is  evident  that  the 
two  peaks  require  almost  three  times  as  many  buses  as  are  needed  during 
the  slowest  part  of  the  day,  and  the  additional  buses  in  service  during 
the  peak  periods  will  have  short  trips,  while  others  may  run  almost  twenty- 
four  hours.  As  far  as  labor  management  agreement  is  concerned  most  of 
them  require  the  same  restrictions  and  basic  demands. 

The  following  is  the  brief  list  of  the  constraints  and  restrictions 
as  founded  in  the  labor  management  agreement  of  one  of  the  transit 
companies: 

The  work  day  of  all  regular  trolley  bus  and  motor  coach  operators 
shall  be  eight  (8)  hours,  in  that  no  regular  run  shall  pay  less  than 
eight  (8)  hours.  Time  and  one-half  shall  be  paid  for  all  work  done 
before  or  after  the  schedule  time  of  the  regular  run  for  all  work 
beyond  eight  (8)  hours  daily.  Forty-six  percent  (46%)  of  all  runs 
shall  be  straight. 

Spread  penalty  shall  be  paid  after  a  spread  time  of  eleven  (11) 
hours  for  regular  operators  and  twelve  and  one-half  (12J0  hours 
for  extra  operations  on  tripper  runs  where  the  intervening  time 
between  taking  out  parts  of  the  run  amounts  is  one  (1)  hour  or 
less,  such  intervening  time  shall  be  paid  as  part  of  run.  Five  (5) 
minutes  pay  shall  be  allowed  operators  making  turn-in  and  will  be 
considered  as  part  of  the  regular  run. 

With  the  above  restrictions,  the  schedule  maker  proceeds  to  develop 
different  possible  runs.  Of  course,  the  best  possible  schedule  for  a  com- 
pany would  be  with  all  straight  runs,  but  due  to  the  variations  in  trip 
lengths,  this  can  never  be  achieved. 


The  use  of  computer  techniques  in  solving  problems  of  scheduling 
men  and  machine  in  public  transit  companies  is  feasible  and  would  be 
simple  were  it  not  for  the  restrictions  which  result  from  the  labor 
management  agreement  and  the  variations  in  the  traffic  situations. 

The  basic  advantages  of  computer  systems  namely  fast  rate  of  pro- 
cessing data  and  accuracy  in  computation  make  the  use  of  computer  pro- 
gramming an  effective  and  better  substitute  for  the  manual  methods 
which  involve  profusive  clerical  work  and  calculations.  The  IBM  1620 
digital  computer  was  used  in  this  thesis.  FORTRAN  (Formula  translation 
system)  which  utilizes  an  automatic  coding  system  resembling  the  language 
of  mathematics  was  used  in  the  development  of  straight  runs. 

Basically,  all  the  program  does  is  simulate  the  motion  of  a  bus 
internally  on  the  computer.  Each  bus  is  followed  across  its  own  route 
and  decisions  and  calculations  are  made  by  the  computer  on  the  basis  of 
the  information  fed  in  advance.  The  same  work  is  dore  in  the  manual  tech- 
nique with  pencil  and  paper  but  would  evidently  take  much  more  time  when 
performed  by  human  beings.  The  restrictions,  such  as  those  imposed  by 
minimum  and  maximum  hours  of  work,  overtime  rates,  spread  penalties, 
and  times  passing  the  relief  point  on  the  route  where  decision  must  be 
made,  are  incorporated  in  the  program  and  implemented  automatically  by 
the  computer.  ; :    , 

For  daily  scheduling  a  two  step  procedure  is  used  which  is  same  in 
both  the  manual  and  the  computer  approaches.  The  steps  are: 

1.  Developing  all  possible  straight  runs. 

2.  Constructing  split  runs  by  combining  pieces  of  work. 

The  main  reason  for  following  a  two-step  approach  is  that  the  computer 
storage  capacity  limits  the  use  of  one  program  to  do  all  the  steps. 
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DEVELOPMENT  OF  STRAIGHT  RUNS 

Step  1; 

The  computer  technique  of  developing  straight  runs  is  the  same  as 
the  manual  technique.  The  headway  for  each  route  supplies  all  the  de- 
tails about  the  input  data.  Such  data  for  Step  1  is  the  route  and  block 
number,  the  pull  out,  pull  in,  and  the  relief  times  for  each  block. 
Cards  are  punched  for  each  block.  The  route  number,  the  block  number, 
and  the  relief  times  are  punched  on  each  card.  If  the  number  of  relief 
times  is  more  than  what  could  be  accommodated  on  one  card,  then  another 
card  would  be  used  for  the  same  block.  All  the  data  received  from  the 
Company  are  in  hours  and  minutes.  The  details  of  Route  26  and  3  Block 
are: 

Pull  out  of  garage  at:  5:25  a.m. 
Pull  into  garage  at:  5:57  p.m. 
Relief  times:  5:25,  7:44,  8:59.  10:40,  11:10 

12:28,  12:58,  2:21,  2:51,  4:12 

4:45,  5:52,  5:57. 

In  order  to  make  the  data  compatible  for  use  on  the  computer  it  had 

to  be  converted  onto  a  24-hour  clock  and  into  hours  and  decimals.  This 

requirement  arises  from  the  fact  that  the  digital  computer  uses  the  digital 

system  of  calculation,  A  conversion  program  was  used  for  this  purpose. 

The  converted  output  for  the  Route  26  and  Block  3  was  punched  out  in  the 

following  format: 

No. /Reliefs  Travel      Route      Block 

13  0.0         26        3 

P/Out  P/In      Reliefs 

5:41  17.95      5:41,  7:73,  8;98,  12:46, 

12:96,  16.75,  17.86,  17.95. 


Once  the  input  data  is  prepared  in  the  manner  prescribed  above,  the 
machine  is  ready  to  start  developing  straight  runs.  The  steps  followed 
in  the  constructing  straight  runs  are: 

1.  The  machine  finds  the  total  trip  time  for  the  block  and  compares 
this  with  eight  hours,  the  minimum  permissible  working  time  for  a 
straight  run.  The  total  trip  time  for  Route  26,  Block  3  is  17.95- 
5.41  =  12.54  hours. 

2.  If  the  total  trip  time  is  more  than  16  hours,  two  straight  runs 
can  be  developed,  one  from  the  front  of  the  block  and  another  from 
the  tail  end  of  the  block,  leaving  a  piece  of  work  in  between.  The 
logic  behind  this  procedure  is  to  avoid  a  late  piece  of  work  which 
would  be  difficult  to  use  in  the  next  step. 

If  the  total  trip  time  is  larger  than  seven  hours  or  less  than  sixteen 
hours,  a  straight  run  and  a  piece  of  work  will  be  developed. 
If  the  total  trip  time  is  less  than  seven  hours,  no  straight  run 
is  developed  but  a  piece  of  work  will  be  punched  out. 
It  was  found  that  for  any  block  having  a  total  trip  time  of  exactly 
seven  hours,  it  would  be  cheaper  for  the  company  to  pay  guarantee 
time  and  make  a  straight  run  rather  than  paying  spread  time  penalty 
and  making  a  split  run.  Because  if  the  company  considers  it  as  a 
piece  of  work  then  in  order  to  make  it  a  split  run  another  piece  of 
work  of  one  hour  will  be  needed.  It  is  a  restriction  in  the  contract 
that  the  minimum  pay  for  a  piece  of  work  has  to  be  at  least  two  hours. 
So,  if  we  combine  this  seven  hours  piece  of  work  in  the  two  hours 
piece  of  work  then  a  payment  of  one-half  hour  of  overtime  has  to 
be  made.  Thus,  total  pay  for  this  run  will  be  9.50  hours.  There 
will  be  a  possibility  that  the  spread  penalty  might  also  be  involved 
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because  of  the  unavailability  of  the  pieces  of  work  within  the 
spread  restriction.  So  it  would  be  cheaper  for  the  company  to  pay 
guarantee  time  and  make  a  straight  run  rather  than  paying  spread 
time  and  overtime  penalties  and  making  a  split  run.  Hence,  it  was 
the  company's  policy  to  use  a  seven-hour  limit  in  deciding  whether 
to  make  a  straight  run  or  a  piece  of  work.  The  total  time  for  the 
example  falls  in  the  second  category. 

3.  The  platform  time  for  the  straight  run  is  now  computed  -  the 
union-management  agreement  states  that  the  company  has  to  consider 
the  five  mintues  turn-in,  travel  time,  and  five  minutes  travel 
allowance,  if  any,  as  a  part  of  the  regular  time.  Therefore,  to 
avoid  paying  overtime,  the  platform  time  is  computed  as  eight  hours 
minum  turn-in  time  (five  minutes),  travel  time  and  travel  allowance, 
if  any.  For  Route  26  and  Block  3,  used  in  the  example,  there  is  no 
travel  time.  Therefore,  the  platform  time  equals  7.91  hours  (8.00 
less  five  minutes  turn-in  time).   The  machine  now  checks  the  relief 
times  on  the  block,  from  the  pull  out  side,  and  checks  for  one  that 
breaks  the  block  into  two  pieces,  one  of  them  being  either  equal  to 
or  slightly  larger  or  slightly  smaller  than  the  platform  time.  In 
the  example,  if  the  block  is  broken  at  relief  12,96  we  get  a  piece 
having  a  platform  time  of  7.55  hours,  where  as  the  next  relief  time 
14.35  gives  a  platform  time  of  8.94  hours. 

4.  The  next  step  is  to  decide  which  of  these  two  relief  times  to 
select  as  the  pull  in  time  of  the  straight  run.  This  selection  is 
chosen  on  the  basis  of  cost.  In  the  case  of  relief  time  12.96,  the 
company  has  to  pay  0,36  hours  as  a  guarantee  time,  but  relief  time 
14.35  pays  an  overtime  of  0.56  hours.  Of  these  two  alternatives. 
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the  one  most  economical  is  chosen.  Therefore,  relief  time  12.96 
becomes  the  pull  in  time  for  the  straight  run. 

5.  The  machine  next  checks  the  trip  time  of  the  remainder  of  the 
block.  In  our  example  this  is  4.99  hours  (17.95  -  12.96  =  4.99 
hours).  It  is  not  enough  for  another  straight  run  so  a  piece  of 
work  is  punched  out. 

6.  Having  developed  a  straight  run  and  a  piece  of  work  from  the 
forward  direction  the  machine  now  constructs  another  straight  run 
from  the  tail  end  of  the  block.  This  is  done  to  give  the  schedule 
maker  the  choice  between  selecting  an  early  or  a  late  straight  run 
on  the  same  block. 

The  straight  run  output  for  Route  26  and  Block  3  is  shown  as 
follows: 
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CONSTRUCTING  SPLIT  RUNS 
Step  2: 

The  next  step  in  constructing  the  schedule  is  that  of  combining 
the  pieces  of  work.  To  explain  this  step,  the  following  two  pieces  of 
work  will  be  used: 

No.   Route   Block    TF    P/Out    P/ln    TB    Platform 

1  1      12     0     6.06    9.51   0       3.45 

2  4       4     0     14.06    18.26   0       4.20 

We  have  seen  above  that  the  pieces  of  work  are  punched  out  one  per 
card,  with  the  following  information:  Route  and  Block  Numbers,  travel 
time  (if  the  piece  begins  away  from  the  garage),  pull  out  time,  pull  in 
time,  travel  time  of  the  pieces  of  work  relieved  on  the  road,  and  the 
platform  time. 

The  computer  program  for  Step  2  is  developed  so  that  the  machine 
performs  the  following  steps: 

1.  The  machine  will  read  each  card  which  contains  all  the  above 
information  about  the  piece  of  work  and  will  store  this  information 
in  its  memory. 

2.  The  computer  selects  the  first  piece  of  work  from  all  the  pieces 
of  work  fed  into  the  machine  and  checks  it  against  all  the  remaining 
pieces  of  work  in  the  same  sequence  in  which  they  are  stored  in  mem- 
ory. This  checking  will  facilitiate  the  development  of  the  split 
runs  by  the  combination  of  two  or  more  of  these  pieces  of  work  within 
the  restrictions  imposed  in  the  computer  program.  These  are: 

a.   A  minimum  gap  of  five  minutes  plus  travel,  if  any,  must  exist 
between  the  pull  in  time  of  piece  (1-12)  and  the  pull  out  time  of 


piece  of  work  (4-4),  if  any.  The  gap  between  the  two  pieces 
being  used  is  4.55;  well  beyond  the  minimum  limit. 

b.  The  spread  limit  is  computed.  There  is  a  maximum  limit 
of  14.0  hours  on  spread  time.  However,  any  time  beyond  11.0 
hours  is  paid  for  at  one  and  a  half  times  the  regular  rate. 
Spread  time  =  (pull  in  time  of  piece  (4-4)  -  Travel  +  turn  in)  - 
(pull  out  time  of  piece  (1-12)  -  travel  +  turn  in)  or  (18.26  - 
0  +  8.3)  -  (6.06  +  0  +  8.3)  =  12.20. 

This  figure  of  12,20  hours  is  within  the  maximum  limit  of  14,0 
hours. 

As  the  above  two  pieces  satisfy  all  the  restrictions,  a  split  run  is 
developed.  The  computer  output  is  in  the  following  format: 


15 


i   .. 


o 

>» 

o 

(0 

• 

a. 

CO 

e 

►-H 

nO 

B 

»— 4 

u 

• 

3 

o 

H 

♦J 

f-H 

o 

C 

^H 

<u 

• 

CL 

o 

s 

^ 

o 

in 

o 

«M 

rr 

(M 

4-> 

• 

• 

CO 

CO 

rr 

1— 1 

a, 

T3 

o 

(0 

rt 

<U 

• 

>-< 

CM 

c 

^H 

C/5 

03 

H 

e 

w^ 

>o 

t-t 

in 

CM 

■N^ 

• 

« 

CL 

o^ 

CO 

1—4 

«J 

3 

-o 

vO 

o 

o 

o 

"^'..^^ 

• 

• 

a. 

vO 

1— « 

tx. 

H 

o 

o 

^ 

u 

o 

CM 

T 

f* 

I— 1 

C3 

0) 

<-» 

3 

»-i 

■^ 

O 

q: 

• 

o 

r-< 

z 

16 

In  a  similar  manner,  the  machine  picks  one  piece  at  a  time  and 
constructs  all  possible  split  runs  combinations  in  the  other  remaining 
pieces.  The  schedule  maker  has  not  to  make  a  selection  from  the  split 
runs  developed. 
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SPLIT  RUN  SELECTION 

Since  the  computer  program  in  previous  sections  is  designed  in 
order  to  construct  all  possible  split  runs,  a  selection  operation  is 
required  by  the  schedule  maker.  Such  a  selection  would  depend  on  the 
objective  of  the  particular  company,  which  sometimes  might  be  governed 
or  restricted  by  the  labor  management.  As  an  example,  some  companies 
aim  at  reducing  the  number  of  trippers  (very  small  pieces  of  work  are 
called  as  trippers),  others  might  have  the  objective  of  combining  the 
large  pieces  of  work  first  and  then  combining  as  many  as  possible  of  the 
left  over  pieces.  Various  methods  have  been  found  for  selecting  the  split 
runs  from  the  output  of  Section  2.  For  this  experiment,  the  computer 
approach  for  the  selection  of  split  runs  has  been  adopted.  The  computer 
approach  of  selection  is  based  on  the  total  pay  time  of  each  split  run. 
The  total  pay  time  of  a  split  run  includes  all  the  allowances,  spread  penalty, 
overtime  penalty,  and  the  trip  time  of  the  constituting  pieces  of  work. 

The  input  data  prepared  for  this  selection  consists  of  cards  punched 
out  for  each  split  run  from  the  previous  section  having  the  information  as 
the  run  number,  route  numbers,  block  numbers,  of  both  the  pieces  of  work 
making  the  run  and  the  total  pay  time  of  the  split  run.  First,  the  machine 
is  instructed  to  arrange  the  cards  in  the  ascending  order  of  the  total  pay 
time  (first  card  for  minimum  pay  time,  second  card  for  little  more  time 
and  subsequently,  they  are  arranged  in  the  increasing  time  order).  Then 
the  machine  picks  up  the  first  card  and  compares  the  route  numbers  and 
block  numbers  of  the  pieces  of  work  for  this  run  with  the  route  numbers 
and  block  numbers  of  pieces  of  work  of  other  runs.  If  the  same  route 
number  and  block  number  are  on  any  of  the  cards  it  will  reject  that  card 
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and  read  the  next  one.  It  thus  goes  on  eliminating  other  runs  and 

finally  makes  a  total  selection  of  split  runs 

in  which  no  piece  of 

.  work  is  used  more  than  once.  Then  it  punches 

out  the  selection  output 

according  to  a  prescribed  format  giving  route 

numbers  and  block  nura- 

bers  and  total  pay  time  for  each  split  run  selected. 

Finally  a  complete  list  of  the  various  split  runs  is  then  given 

when  the  minimum  total  pay  time  is  selected. 

■■'-  ■  ';  •  '          ■■  -,  ''''".      •^'''''"■•■'  --  ^ 
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SELECTION  OF  PIECES  OF  WORK  NOT 
USED  FOR  SPLIT  RUNS 

After  making  the  selection  of  split  runs  it  becomes  necessary 
to  eliminate  the  pieces  of  work  used  for  making  the  split  runs  from 
the  total  pieces  of  work  used  for  Section  2.  The  main  idea  behind 
this  process  of  elimination  of  pieces  of  work  used  is  to  allocate  the 
penalty  for  the  pieces  which  could  not  be  used  in  Section  2.  This  is 
done  in  the  following  steps  by  making  use  of  the  computer. 

(1)  The  output  of  "Split  Run  Selection,"  is  stored  in  the  machine. 

(2)  Then  all  the  pieces  of  work  from  Section  2  are  read  by  the 
machine  one  by  one,  and  each  piece  of  work  read  by  the  machine  is 
compared  with  the  pieces  of  work  which  are  used  for  selected  split 
runs  stored  in  the  machine's  memory.  If  the  piece  of  work  read 

by  the  machine  is  not  used  for  the  split  runs,  then  it  is  punched 
out.  Otherwise,  it  reads  another  piece  of  work  and  thus  again 
compares  it  with  the  pieces  of  work  used  for  split  runs. 

(3)  The  Steps  1  and  2  are  repeated  until  all  the  pieces  of  work 
are  read  and  compared  with  pieces  of  work  of  split  runs. 

The  output  is  punched  in  the  same  format  as  of  the  piece  of  work 
of  Section  2.  ^S '^ 
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ALLOCATING  PENALTIES  TO  THE  UNUSED 
PIECES  OF  IVORK 

The  pieces  of  work  which  are  not  used  for  the  development  of 
straight  runs  are  then  allotted  a  penalty  for  the  purpose  of  calculating 
total  pay  time  of  the  company  for  a  particular  division.  The  system  of 
assigning  a  penalty  for  the  pieces  of  work  not  used  differs  from  one 
company  to  another.  An  usual  system  of  assigning  a  penalty  adopted 
for  carrying  out  this  research  is  given  below: 

PLATFORM  TIME  OF  PIECE  OF  WORK  PENALTY 

BETIVEEN   190.  Min  -  299  Min  4  hours 

300.  Min  -  499  Min.  ^6  hours 

.      500,  Min  -  Above  8  hours 

By  adopting  the  above  penalty  system,  the  total  time  for  the  pieces 
of  work  not  used  is  then  calculated. 

Generally,  the  trippers  (pieces  of  work  having  platform  time  less 
than  190.  Min)  are  not  assigned  any  penalty  by  the  transit  companies.  So 
tripper  times  are  just  added  together  which  will  give  total  tripper  time 
for  the  complete  division. 

The  total  time  for  which  the  company  makes  payment  to  the  operators 
daily  will  consist  of: 

1.  Straight  run  pay  time. 

2.  Split  run  pay  time  (including  spread  and  overtime  penalties). 

3.  Penalty  for  pieces  of  work  not  used  for  split  runs. 

4.  Total  tripper  time. 


^,r'  J-  ^mr^*-. 


22 
DESIGN  OF  EXPERIMENT 

Here  we  are  concerned  with  the  establishment  of  the  mathematical 
or  statistical  relationship  existing  between  a  number  of  economic  variables. 
A  model  or  hypothesis  as  to  the  assumed  relationships  between  the 
economic  variables  is  to  be  constructed.  Then,  economic  measurements 
are  applied  to  each  variable  and  the  degree  of  relationship  is  determined. 
There  are  various  constraints  and  restrictions  in  labor  management 
agreements  of  public  transit  companies.  There  are  always  a  few  constraints 
like  platform  time  and  spread  time  in  the  agreement  on  which  controversy 
is  always  going  on  between  the  trade  union  and  the  management.  While 
conducting  this  experiment  for  predicting  total  pay  time  of  operators 
in  public  transit  companies  it  was  felt  that  platform  and  spread  time  are 
the  two  main  deterministic  variables.  Here  we  are  establishing  the  likeli- 
hood that  these  two  variables  have  a  relationship  to  the  total  pay  time. 
So,  the  model  building  for  the  prediction  of  total  pay  time  seeks  the 
basic  pay  time  determinants  such  as  platform  time  and  spread  time. 

In  Step  1  of  straight  run  development  there  is  a  restriction  that 
time  and  one-half  is  paid  as  overtime  for  all  work  beyond  eight  hours. 
Also,  in  Step  2  of  split  runs  development  there  is  a  restriction  that 
the  spread  time  must  not  exceed  fourteen  hours.  So  the  two  constraints 
which  are  varied  for  this  experiment  are  spread  time  and  platform  time. 
The  platform  times  are  varied  from  700  minutes  to  900  minutes  and  the 
spread  time  is  varied  from  1150  raintues  to  1450  minutes.  The  input  data 
used  for  the  various  sets  of  combination  of  spread  time  and  platform  time 
is  obtained  from  one  of  the  public  transit  companies.  The  sequence  of 
steps  is  as  follows: 


-"TS^^nrrrrTTr^ 
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(1)  Development  of  straight  runs  as  done  in  Step  1  of  Scheduling, 
introducing  various  platform  time  limits  say  700  minutes,  725  minutes, 
850  minutes  and  on. 

(2)  Selection  of  straight  runs,  trippers  (pieces  having  platform 
time  less  than  190  minutes)  and  pieces  of  work. 

(3)  Development  of  split  runs  from  the  pieces  of  work  as  done  in 
Step  2  of  Scheduling  by  introducing  various  spread  time  limits  say  1150 
minutes,  1200  minutes,  1250  minutes,  and  so  on, 

(4)  Selection  of  split  runs  on  the  basis  of  total  pay  time  which 
mainly  includes  spread  penalty,  overtime  penalty  and  the  trip  time  of 
pieces  of  work.        ' 

(5)  Making  a  list  of  the  pieces  of  work  not  used  for  the  development 
of  split  runs  and  assigning  penalties  to  them. 

(6)  Making  a  list  of  trippers  and  adding  their  trip  times. 

Figure  4  gives  the  tabulated  results  for  total  pay  time  in  detail 
for  the  various  platform  times  and  spread  times.  The  total  pay  time  for 
operators  of  a  complete  division  consists  of  time  of  straight  runs,  time 
of  split  runs,  penalty  time  for  pieces  not  used  for  split  runs  and  total 
tripper  time. 
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MXTIVARIATE  AND  BIRARIATE 
CORRELATIONAL  ANALYSIS 

While  conducting  an  experiment  for  predicting  total  pay  time  of 
operators  in  public  transit  companies  it  is  felt  that  platform  time 
and  the  spread  time  are  the  two  deterministic  variables.  Here  we  are 
establishing  the  likelihood  that  these  two  variables  have  a  relation- 
ship to  the  total  pay  time.  The  hypothesis  that  Y  =  f(Xj,  Xg),  namely 
that  Y  is  a  function  of  x,  and  Xg  where, 

Y  =  total  pay  time 

x,=  platform  time 

X2=  spread  time. 

The  data  presented  in  Figure  4  for  the  total  pay  time  (Y)  together 
with  two  dependent  variables,  platform  time  (x. )  and  spread  time  (X2)  is 
used  for  establishing  correlational  equations.  Both  spread  time  and 
platform  time  are  expressed  in  minutes  and  fraction  of  minutes.  The  total 
pay  time  is  expressed  in  hours  and  fraction  of  hours. 

A  Symbolic  Programming  System  program  has  been  used  for  finding  the 

regression  equation.  The  equation  is  developed  in  the  multiplicative 

2 
form  and  the  coefficient  of  multiple  determination  (R  )  is  also  deter- 
mined which  enables  one  to  know  how  close  the  estimated  values  are  to 

the  actual  values. 

2 
The  standard  formula  for  R  ,  taken  from  Nemmers  (1962)  is. 


^y       2  V     o-J'-N-M-^ 
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)  =  The  sum  of. 

d  =  The  deviation  of  the  actual  values  from  the  estimated. 
N  =  Number  of  observations. 
Y  =  The  actual  values. 

M  =  The  number  of  constants  in  the  multiple  regression  equation. 
The  first  predicting  equation  is: 

log^  (Y)  =  9.58359  +  0.05147  log^  (x.)  +  0.41229  log  (x^) 


or 


V  _  9,58359  ,     ,0.05147  ,     .-0.41229  ,,, 

Y  -  e      .  (x.)      ,  (xg)  (1) 


where   x.  =  Platform  time  in  minutes. 
X2  =  Spread  time  in  minutes. 
Y  =  Total  pay  time  in  hours  and  fraction  of  hours. 

2 

For  the  above  equation  R  is  0.87659  or  R  =  0.94  which  shows  that 

there  is  not  much  explained  and  unexplained  variation.  By  putting  the 
different  values  of  x.  and  x^  in  the  above  equation  Y's  are  computed  which 
are  tabulated  in  Figure  5,  and  graphs  are  plotted  for  Y  against  x,  and  Y 
against  Xg,  which  shows  that  the  total  pay  time  decreases  by  about  10%  for 
the  spread  time  change  from  1125  minutes  to  1450  minutes.  Also,  it  is 
seen  that  for  change  in  platform  time  from  700  minutes  to  900  minutes, 
the  total  pay  time  increases  by  about  1*5%. 
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The  previous  equation  for  the  total  pay  times  expressed  as  a 
function  of  spread  time  and  platform  time  is  in  the  multiplicative  form 
and  it  gave  a  coefficient  of  multiplicative  determination  of  0.877,  Then 

an  attempt  was  made  to  change  the  form  of  equation  expressing  it  as  a 

2 

function  of  the  same  variables  so  that  the  value  of  R  may  become  greater 

than  the  previous  value.  The  following  is  the  second  form  of  the 
predicting  equation  for  the  total  pay  time  in  terms  of  the  two  deter- 
minants platform  time  and  spread  time. 

Y  =  4207.75112  -  445-39419  log^  (X2)  +  0.06529  (x^)  (2) 

where  Y  =  Total  pay  time  in  hours  and  fraction  of  hours. 
x,=  Platform  time  in  minutes. 

X2=  Spread  time  in  minutes. 

For  the  above  equation  (2)  the  coefficient  of  multiplicative  deter- 

2 
mination  (R  )  is  0.903  or  R  =  0.95.  It  can  thus  be  concluded  that  a 

better  correlation  of  spread,  plat  and  total  pay  time  is  given  by  this 

equation.  From  this  above  equation,  Y's  are  calculated  for  different 

values  of  x,  and  X2  and  then  they  are  tabulated  in  Figure  6.  The  graphs 

are  also  plotted  for  Y  against  x.  and  Y  against  X2.  It  is  seen  that 

total  pay  decreases  by  about  9*8%  for  change  in  spread  time  from  1125 

minutes  to  1450  minutes.  Also,  total  pay  increases  by  about  12%  for  a 

change  in  platform  time  from  700  minutes  to  900  minutes. 

On  analyzing  both  the  equations  which  give  total  pay  time  as  a 

function  of  platform  time  and  spread  time  it  is  noticed  that  total  pay 

time  is  not  as  significantly  affected  by  platform  time  as  it  is  by  spread 
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time.  So  for  each  platform  time  the  total  pay  time  is  expressed  as  a 
linear  function  of  spread  time  and  the  coefficient  of  multiplicative 
determination  is  determined.  The  results  are  summarized  as  follows: 


Platform  Time 
(Minutes) 

Equation  Y  =  Total  Pay  Time  (Hr.s) 
X2=  Spread  Time  (Mins.) 

8^ 

700 

Y  =  1519  -  0.346  X2 

0.983000 

750 

Y  =  1476  -  0.331  X2 

0.997030 

800 

Y  =  1458  -  0.293  x^ 

0.996660 

825 

Y  =  1517  -  0.349  X2 

0.995060 

850 

Y  =  1580  -  0.383  X2 

0.999631 

875 

Y  =  1533  -  0.358  X2 

0.999495 

900 

Y  =  1473  -  0.306  X2 

0.995820 

To. confirm  statistically  that  platform  time  does  not  effect  total 
pay  time  as  significantly  as  spread  time,  two-way  analysis  of  variance 
is  tested  out  in  the  next  section. 
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TWO-WAY  ANALYSIS  OF  VARIANCE 

What  is  here  to  investigate  is  the  effect  of  K,  the  different 
platform  time  restrictions,  and  n,  the  different  spread  restrictions 
on  the  total  pay  time.  Using  notation  introduced  in  Freund  (1963),  let 

the  "Plats"  be  A.  for  i  =  1,  2 k  the  spreads  be  B.  for  j  =  1,  2,.,.n. 

Suppose  that  total  pay  time  is  found  for  each  possible  combination  of 
plats  and  spread,  that  the  (n,  k)  different  pay  times  are  randomized, 
and  that  x. .  is  the  pay  time  obtained  with  the  i   plat  and  j   spread. 
A  possible  model  for  this  kind  of  problem  is  to  look  upon  the  x. .  as 

values  assumed  by  independent  random  variables  having  normal  distributions 

2 
with  the  means  u. .  and  variance  a  ,   where 


and 


B 


2^   a.=0   and  2]  ?■  =  0 
i=l  .   j=l 


where  a.  is  the  effect  of  i   plat  and  p.  is  the  effect  of  j   we  could 
also  specify  these  assumptions  by  writing  spread, 

x..=^+  a.  +  p  +  e   for  !:}'^"--|;  (4) 

ij        1    J    ij  .   j-l,<i,...n 

where  e.  .  are  values  assumed  by  independent  random  variables  having  normal 

2 
distributions  with  0  means  and  the  common  variance  a  ,     It  is  to  be  seen 

that  in  this  model  the  effects  of  the  two  variables,  that  is,  the  a.  and 

•       •      1 

p. ,  are  added  to  n. 

The  null  hypotheses  we  shall  be  testing  are: 
(1)  That  the  a.  are  all  equal  to  0,  and 
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(2)  That  the  p.  are  all  equal  to  0,  the  corresponding  alternative 

hypotheses  are  that  the  respective  parameters  are  not  all  equal  to  0. 

The  tests  of  these  hypotheses  are  based  on  the  following  analysis  of 

the  total  variability  of  the  data,  decomposing  it  into  terms  attributed 

to  differences  among  the  A's,  differences  among  the  B's,  and  chance 

(experimental  error):  . 

k       n  k  •     n 

yy         -2     y--2y..2 

L     L      (x. .  -  x)    =  n    Zj        (x.  .  -  x)     +  k  Z_i      (x.  .  -  x) 
i=l  j=l        *J  i=l  ^  j=l  J 

k       n  .  (5) 

Y      Y  -  -  -  2 

L     L      (x. .  -  x..  -  x, .  +  x) 

i=l  j=l         iJ         1   .         J 
where 

X. .  =  -     (  )      X.  .) 

i=l  •   '    ■ 

n 

X..  =  i  (     Y     X..) 

in        Zj      ij       . 

j=l  ■    ■     •    '■        - 

k       n  •  ■ 

X     =  -r-   (Y    r  X..) 

k.n        iLi      Li.      ij 


-.j^yj- 


Equation  (5)  may  be  written  as, 

SST  =  SSA  +  SSB  +  SSE  ; 

where 

SST  =  Total  Sum  of  squares. 

SSA  =  Sample  sum  of  squares  for  variable  A. 

SSB  =  Sample  sum  of  squares  for  variable  B, 

SSE  =  Error  sum  of  squares. 
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It  is  seen  that  if  the  null  hypothesis  concerning  the  o.  is  true, 

2        2 
then  SSA/a  and  SSE/a  are  values  assumed  by  independent  random  variables 

having  chi-square  distributions  with  k-1  and  (n-1)  and  (k-1)  degrees  of 

freedom;  if  this  null  hypothesis  is  not  true,  then  SSA  can  be  attributed, 

at  least  in  part,  to  differences  among  the  A's,  that  is  differences  among 

the  platform  time.  Similarly,  if  the  null  hypothesis  concerning  the 

2         2 
p.  is  true,  it  can  be  seen  that  SSB/a  and  SSE/a  are  values  assumed  by 

independent  random  variables  having  chi-square  distributions  with  (n-1) 

and  (n-l)(k-l)  degrees  of  freedom;  if  this  null  hypothesis  is  not  true, 

then  SSB  can  be  attributed,  at  least  in  part,  to  differences  among  the 

B's,  that  is,  differences  among  the  spread.  If  both  of  the  null  hypothesis 

are  true,  it  is  seen,  furthermore,  that  SST/a  is  a  value  assumed  by  a 

random  variable  having  a  chi-square  distribution  with  (nk-1)  degrees  of 

freedom: 

In  accordance  with  the  following  theorem: 

"if  X.  and  X2  are  independent  random  variables  having  chi-square 
distributions  with  v,  and  Vg  degrees  of  freedom,  then; 

■    '- v^  •    • 

V^2     •■•■,■.::  *  ■  ;,  ■        -^, 

has  an  F  distribution  with  v.  and  v^   degrees  of  freedom. 

The  test  of  the  null  hypothesis  concerning  the  a.  can  be  based  on 

the  static: 

p  -  SSA/(k-l)      ,  (n-1) (SSA)  ... 

■^A  ~  SSE/(n-l)(k-l)  ■     SSE       *•*•  ^^^ 

which,  under  the  null  hypothesis  that  the  a.  and  all  equal  to  0,  is 

a  value  assumed  by  a  random  variable  having  the  F  distribution  with 

k-1  and  (n-l)(k-l)  degrees  of  freedom.  We  reject  this  null  hypothesis 
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if  F.  is  greater  than  or  equal  to  the  critical  value  given  in  Table  VI  b 
of  Freund  (1963).  -, 

Similarly,  the  test  of  the  null  hypothesis  that  the  p.  are  all  equal 

•I 

to  0  can  be  based  on  the  static, 

P  _  p.  _  SSB/(n-l)      _  (k-1)  SSB 
''b  ~  ''b  ~  SSE/(n-l)(k-l)       SSE 


(7) 


which,  under  the  null  hypothesis  that  the  p.  are  all  to  0,  is  a  value 
assumed  by  a  random  variable  having  the  F  distribution  with  n-1  and 
(n-l)(k-l)  degrees  of  freedom.  We  reject  this  null  hypothesis  if  Fj,  is 
greater  than  or  equal  to  the  critical  value  given  in  Table  VI  b  of  Freund 
(1963).. 

The  analysis  of  variance  table  for  this  kind  of  a  two-way  analysis 
is  usually  presented  in  the  following  fashion: 


Source  of 
Variation 

Degrees  of 
Freedom 

Sum  of 
Squares 

Mean 
Square 

F 

Between 
A's 

k-1 

SSA 

,,„,_SSA 
"S^k-1 

MSA 
MSE 

Between 
B's 

n-1 

SSB 

MSB=SS? 
n-1 

MSB 
MSE 

Error 

(n-l)(k-l) 

SSE 

MT-    SSE 

"^-  (n-l)(k. 

■1) 

1 

Total 

nk-1 

SSI 

FIGURE  7 
ANALYSIS  OF  VARIANCE  TABLE 


Here  the  mean  squares  are  again  the  sums  of  squares  divided  by  the 
respective  degrees  of  freedom.  •    . 
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The  experiment  designed  to  test  seven  different  spreads  and  five 
different  plats  yielded  the  results  of  total  pay  time  shown  in  the 
following  table:  (Total  plat  is  in  hours.) 


^s.  Spread 
PlatX^^ 

1150 

1200 

1300 

1400 

1425 

5 

&ij 

700 

1127 

1102 

1070 

1022 

1031 

5352 

750 

1100 

1076 

1047 

1003 

1006 

5232 

800 

1128 

1114 

1074 

1047 

1037 

5400 

825 

1121 

1093 

1059 

1025 

1015 

5313 

850 

1140 

1122 

1083 

1032 

1034 

5417 

875 

1120 

1120 

1059 

1022 

1015 

5336 

7 

7x.. 
^  ij 

7863 

7725 

7468 

7194 

7159 

C 

37425  ) 

J 

5  7 


2 


X.  . 

ij 


j=li=l 


FIGURE  8 
ACTUAL  TOTAL  PAY  TIME 


Using  the  modified  form  of  formula  (5)  we  get, 

k       n      •  k       n  2 

Li    'L    ij   "  k*n    LZj      L      ij  J 
i=l  j=l  i=l  j=l 


7       5 


-J   Ly  ■   Trfe,   [  I    Z  ^j  ] 


1=1      j=l 


7       5  2 

X.  . 


i=l  j=l 


1 


=  4066645  -  ^  .    (1400630625) 
=  4066645  -40018018 
=  48627 
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k    n     2 
>..  . 


i=l  JtI  i=l  j=l 


SSB  = 


7    5     2      7   5     2 
i=l    j=l         i=l  j=l 


1  (200113107)   i  (1400630625) 
"5  "35 


=  40022621  -  40018018 
=  4603. 

j=l  i=l  1=1  j=l 


and  SSE 


5    7     2         "^       ^  ', 

7  I  [  I  ^ij  ]   -  fe  [  I  I  '^ij  ] 
j=l  1=1  1=1  j=l 

40039762  -  40018018 

21744 

SST  -  SSA  -  SSB 

48627  -  4603  -  21744 

22280. 


Source  of 
Variation 

Degrees  of 
Freedom 

Sum  of 
Squares 

Mean 
Square 

F 

Between 

A's 

6 

4603 

767.1 

0.827 

Between 
B's 

4 

21744 

5436 

5.85 

Error 

24 

22280 

928 

Total 

34 

48267 

FIGURE  9 
COMPUTATION  FOR  ANALYSIS  OF  VARIANCE 
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Critical  values  of  the  static  F  for  Table  VI  b  of  Freund  (1963) 


are: 

^•01,6.24  =  3.67 

^•0..4.24  =  4.22 

Since  F.  =  0*827  which  is  less  than  F.,,,  ,   ^.   and  F„  =  5'85,  much 
A  01,6,24     B 

greater  than  F.„-  .  ^.,   the  null  hypothesis  for  A's  cannot  be  rejected, 
but  the  null  hypothesis  for  the  B's  is  rejected  at  the  .01  level  of  sig- 
nificance. We  may  conclude  that  the  total  pay  time  is  not  much  affected 
by  platform  restrictions  and  is  significantly  affected  by  different  spread 
restrictions. 


■■^- 
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GRAPH  2 
BETWEEN  TOTAL  PAY  FROM  EQUATION  1  AND  PLAT 
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GRAPH  3 
BETIVEEN  TOTAL  PAY  FROM  EQUATION  1  AND  SPREAD  FOR  VARIOUS  PLATS 
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PRCGRA".  MC.  1      ,   ,   .  • 

CCNVERSIGN  PROGRAM  FOR  HOURS  MINUTES  TO  HOURS  DECIMALS 
DINfiNSIOiM  RnL(lOO) 

Ull  RFADlli2»TRAVEL      '  ■ 
1112  FORMAT  (F?-.0)  ■   ' 

TRAVFL~TIME(TRAVEL»1  .667) 
1 1  1  N  ]  =  1 

r;?=i3 

REA01»LI,\E,N0TRIP»P0UT,PIN,  (REL(  I  )  »I=N1»N2) 
1  FORMAT(2I2»15F5.0) 

IF(LINE) 100,100,2 
7     IF(REL(N2) )3,4,3 
„.   ?.  IF(REL(N2)-PI.\)5»6o 

5  Kl=i+N2 
N2=N2+15 

READ1,LINE,N0TRIP» (REL( I ) > I=Ni»N2) 
GO  TO  2 

4  DO  7  1=1, N2 

IF(REL.(I  )  )7,0,7       .,  '. 
7  CONTINUE  .:...,  ^  :v-  .. 

8  1  =  1-1  '  :     -.„...    ivC;:/ 

GO  TO  9       ■.'.''    ■"' 

6  I=N2  ,  .-  .  ., 
9  DO  10  J=1,I 

m  RFL(J)=TIME(REL( J) ,1.667) 
POUTC= TIME (POUT, 1.66 7) 
POINC=TIME(PIN,  1.667) 

17  PUNCH18,  I  ,TRAVEL»LIN£»NGTRlP»POUTC,POIi\C,(RbL(  J)  ,J=1,I  ) 

GO  TO  111  ''-.v  ::  ■  ' 

lUO  PAUSE  ■    ;    ■:...- 

GO  TO  1111       '      •■ 

18  FORMAT ( I3,F5.1,2I3,llF6.0/26X,9F6.0/26X,9F6.0/26X,9F6.n/26X,9F6.0/ 
1?6",,9F6.C'/26X,9F6.0) 

END  '   '  • 

*■»!•  TIME  FUNCTION  ■  .-  - 

FUNCTION  TIME(X,Y) 
IF(X)1 ,2,2 

1  X  =  -X+1200.   '-'■".'■' 

2  ix  =  x/loo.        :  ■.: 

HR=IX*10U 

TlME=(X-HR)«Y+HR  •- 
RETURN 

END     .  ;    : 


wTTrygigt'jviJtBeiqe 
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PROGRAM   NO.  2.  " 

STRAIGHT   RUN  MAKING  PROGRAM 


DIMENSION   RELdOC) 
2  TURN=8.33 


TR  =  0 
2'-^].  RFL{  n=C 
10  DO  201  Ii=l  ,100 

PU^'CH  20 2 
2^2  FORMAT (//  )    '  ■•  ■  ■ 

R  EAD 11  ,  N'i ,  T  RAV  EL  ,  ROUTE  ,  SLOCK  ,  T  I  MOi\  ,  11 MOFF  ,  {  R  EL  (  L  )  .  L=  1  ,  M  ) 
lA  aACK  =  C      .      ..r,^._. 

FRONT  =  '.-   ■     .    "-       ■  ■  '■    ■  ■ 

POUT  =  Tlf-''ON 

POIM  =  TP-',OFF 
IB    FIGHT=600.-TURN 
?"■    A I  LOV'=0 

21  IF(TIK0FF-TIM0N-6C0. )?4,24,22 

22  B/\CK  =  1 

2  4  I  F  {  T  I  MOFF-T  I  tV;ON  )  1 0  , 1 0  ,  2  5 
25  TRIPT  =  Tiy.OFF-TIMON' 
2  6  IF(TIM0N-P0UT)2  7,2  8,2  7   ' 
27  TF=TRAVEL 

2  8  IF(TIN'OFF-POIK)29,30»29 

29  TR=TRAVFL 

30  IF(TRIPT-6C0. )31»136 »136 
136  IF(TRIPT-5r;0.  )35,36,150  . 
I'^r,     IF(FRONT)90';r.,50,2^0 

31  IF(TRIPT-190. )331>331^43T■ 
33  1  PUNCH  33  2  , ROUTE* SLOCK  ,TF  ,  T  I  MOiN  ?  T  I  MOFF  ,  75  V  TR  I  PI 
3  3  2  F  ORMAT  (  5X  ,  2  I  4  ,  F 5  .  1  »  2  1 6  ,  F  5  .  1  ,  !i X  ,  I  6  ,  29X  ,  i H-  ) 

GO  TO  33- 
431  PUNCH32  ,  ROUTE  ,  BLOCK  ,  TF  ,  T  IMOiX  ,  T  I  MOFF  ,  Tb  ,  TR  IPT 

3  3  TF=C;  ,.  ,  ,  ,  - 

TB  =  C  "■■"  ■;   ,     . 

I F( bACK ) 9000 ,10,100 

36  IF(TIV.O.M-POUT)37,3  8,37 

37  7F=TRAVEL 

3  8  IF(TP-'0FF-P0IN)39,4C>39 

30  tr=travfl    ■      .  • 
40  allow=tf  +  tf<. 

TOTAL=TRIPT+TURM+TF+TB 
42  PU'NCHA  3  »  K,  ROUT  E  » BLOCK,  T  F,  T  1M0N,T  I  MOFF,  To,  TR  IPT,  ALLOW,  TURN*  TOTAL 
TF  =  0 

TURM2=16.67 
K=K+1  :  . 

TB=c  ' ; .. 

I  F{  BACK)  100,10,100  ■  ..   ■- .    . 

FR0NT  =  1  V  -  .:■  . 

5]  DO  5  3  L  =  l  ,'•' 
50  PLAT  =  T  I  M0'N+60r:. -TRAVEL-TURN 
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b? 

1 F(PLAT-REL(L)  ) 65 » 34,  S3 

5 '3 

CCnTINUL 

IF (L-V) 53,3 1.31 

"^A 

TP=TRAV-L 

5  5 

TCTAL  =  PFL(L)-TIMC^'^-TUR^I  +  TRAVFL 

•^6 

AL  1-CW=TF  +  Tf3 

tript=r[:l(l)-tin'on       :- 

i 

K=^:<  +  1                    '  ' 

59 

PUi\CH43  .  K  ,  ROUTE  ,  BLOCK  ,  TF  ,  T  I  MCN  ,  REL  (  L  )  ,  Ta  ,  TR  1  PI ,  ALLOW  ,  TURN  .TOTAL 

FR0NT=1  •   .        .  : 

TF=o          ■■   ; 

TP  =  ^J       ■  ■    ■  : 

TiMOf!  =  REL{L)   ' 

6  3 

EIGHT=600.-TURNtTRAVEL 

rO  TO  2^ 

65 

P'."C=(RFL{L)-PLA'T)*.5 

65 

PK'P  =  PLAT-RFL(L-1  ) 

67 

I F(PWO-PWP) 6  8.6  8.74 

68 

;<  =  K  +  i 

TP  =  TRAVEL      '               ,.  . 

ALL0W  =  TF  +  Tr3 

T0TAL=REL(L)-TIM0N+TURN+TF+TB 

TRIPT  =  REL(L)-TIMO.\ 

GO  TO  59 

74 

K  =  K  +  1                --:•, 

Tp.  =  TRAVFL 

ALL0v;=TF  +  T3+PWP                                                  .    ' 

TCTAL  =  RFL(L-1)-TI^10N  +  TURM+ALL0V/ 

TRIPT  =  REL[L-l)-TIM0^4 

PUNCH43.  <,  ROUTE,  BLOCK.  TF.T  I  MO.M,  RE  L(L-l)  ,  T  B,  TR  I  PT  ,  ALLOW  ,  TURN  ,  TOTAL 

TIM0N=RFL(L-1 > 

TB=l 

FP0NT=1 

GO  TO  63     .  •  .. 

99 

TE=TRAVEL 

100 

Tir-'ON  =  ?OUT  .  • 

TIMOFF  =  POIN! 

20  0 

PLAT  =  TlV:0FF-600.+TRAVEL  +  TURN 

301 

DO  103  L=1.M 

IF{PLAT-RFL(L) ) 121 ,104,103 

2  O'a 

rOMTIMUF 

IF(L-iV)  103,31,31 

K4 

TF=TRAVFL 

ALLOW=TF+TB         "   •   :' 

TOTAL=6C0. 

K.  =  K  +  1 

TRIPT=TIMOFF-REL(L) 

109 

PUNCri  43.K. ROUTE. BLOCK, TF. REL  (L)  ,  T  IMOFF  ,  Tl3,  TR  I  PT  ,  ALLO.V ,  TURi'M  ,  TOTAL 

FRONT  =  o             ...  .  . 

111 

TF=(; 

113 

T!MOFF  =  REL(L)       \      ,.-;  .■ 

t       "'  -         * 

50 


114  TIvi^  =  TlvCFF-TIN'CN 

115  !F(TIf'£)9'-:'0,in,116 

11A     !F(T!'-'F-(6rr>..-TRAVFL-TRAVHL)  )117,]  19,109 

]nQ       T''.  =  TRA\/FL 

OC    Tr;    5^'       . 
1, 1  7    SACK  =  0 

GC    TC    2  5 

T  1  o    n  ArK='^  "  .   .      .  ■ 

120    GO    TO    36 

121  pwc=(plat-u£l{l-1) )-.5 

122  pwp  =  re:l(L)-plat 

12  3  if(pwc-pw?)124,12a->134 
124  tf=travfl 

ALLCW=TF+Tn 

tc;tal=timcff-rel(l-i  )+turn+ allow 

K  =  <  + 1 

TPIPT=TIMCFF-RRL(L-1 ) 

PUNCH  4'3  ♦<»  ROUTE.  BLCCK,TF,REL(L-1  )  ,  T  I  ,V0FF  ,  Tb  ,  TR  I  PT  ,  ALLOW  ,  TURN  »  TCTA 

1 L  ■-•..,• 

FROi\T  =  t^  ■ 

TF=C;  '.  . 

TIM0FF=REL{L-1 ) 
GO  TO  114 
134  TF=TRAV£L 

ALLOW  =  TF+TR+PV.'P 

TOTAL  =  TIMOFF-REL(L)+TUR^!+ALL0W 
•<  =  K  +  1 

TRlPT=TIMOFF-RtL(L) 
GO  TO  109 
oOOO  STOP  "'•'^•-  '  ' 

3  2   FORN^AT  (9X,2l4»F5.1  »2l6.F5.1  .5X,I6) 

4  3  FORMAT (  I  3 > 6X ♦ 2  I 4,F5 . 1 » 2 16 » Fb . 1 > 5X , i 6 »2 F6 • i , F 7 . 0 ! 

11  FCi'<iMAT(l3,F3.1,2F3.0»llF6.0/26X»9H6.0/26X,9F6.0/26X,9r6.0/26X,9F6» 
10/26X,9F6.0)    .  ■•  .  ■, 
END         '■  ■  ':   •  '  ■       •  . 
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PRCCiRAM   MC.  3 

A  SPLIT  RL-'N  DFVFLCPN^'INT  PROGRAM 

DI  .^^t  MSI  CM  ROUTE  (  200  )  ,6LCC;<  (  200  )  »T  F  (  2uO  )  tT  I  N'ON  (  200  )  ,  T  IMCFF  (  200  )  ,  TtJ  ( 

1200)  ,TKIPT(20'.-)           

1  r?EAD2.IP'A 

2  rCRy.AT(3XI3)        '   ■ 

|V;  —  U                                              .                             >-■■-,:■ 

3  Tu;^i\  =  S.33                   ■   ' 

A  FI VE=8.33 

r  ■  =  ^ 
6  DC  9  1  =  1  ,IPW 

o  rO^'TTNUF 

in   L=l               ■     '               .         \ 

n   00  loov,/  j  =  i , ipv/  ■                        •                                                     „     , 

7   Rr AD8 , ROUTE ( I  )  , BLOCK ( I)  ,  TF ( I)  , T I MON ( I ) ,TIMOFF {  I)  ,  T3 ( I)  ,  TR I PT {  I) 

12   PE.VLT1=0   . 

PFNLT2=0 

PENLT3=0 

16  I  F  {  T  I  MON  (  J  )  -  T  I  ••'ON  (  L  )  )  1 000  .  1  000  »  1 7 

17  IF(TIMOFF{  J)-TliVOFF(L)  )1000»lCOO,ia 

18  IF{  TIf^'OiM(  J  )-TIMOFF(Ll  11000,1000  >19 

TO  GAP  =  Tr.'.OM(  J)-TIMOFF(L) 

20  (i'5"A'<:  =  TR(L  )  +  TF(  J)+FI  VE 

21  I  P(r:.AP-o(5c-AK  )  3  00^,??  ,2? 

2  2  I  F(GA  P- 1  O.}.  )  50^,23,2  3 

2  3  SPRFAn)=Tiy'OFF{  J)-TIVOM(L  ) +  TF  (  L  ) +Tt3  ( J  ) +TURN 

?.h    IF(1400. -SPREAD)  1000, 25. 12A 

124  IF(SPREA0-110:'.  )26,26,25 

25  PFNLT3=( SPREAD-1100. ) 

2  5  P Wl  =  TR I PT ( L ) +TF ( L ) +Tb ( L  1 

27  PW2  =  TRIPT(  J)+TF(J)+Tti(J) 

28  IF{PW1-1UG.) 29,31,31 

2  9  PENLT1  =  100.-Pva 

3  0  PW]=10^.                ;'.■,■ 

31  IF(PW2-}or'.)32,34,34 

^2  PFNLT2=100.-PW2 

-J  3  pi,,?  =  10(,. 

34  TCTAL  =  PW1  +  PW?  +  2.*T1)RN 

35  IF(TOTAL-75C. )1000,36,45                                    • 

,  36  PFMLTY=50. 

37  PAYTIM=TOTAL+PENLTY 

ALL0W=PtMLTY+PENLTl  +  PENLT2+  TF(L)  +  io(J)  +TLi{L)+TF{J) 

M=M  +  1 

40  PUMCH41,R0UTE(L)  ,i3L0C.<(L)  ,  TF(L)  ,TIMO,N(L)  ,TIMOFF  (L)  ,Ty(L)  ,TRIPT(L) 

42  PUNCH43,M,R0UTE(J) .BLOCK! J) ,TF( J)  ,TIMON(J)  ,TIMOFF(J) ,TR( J) ,PENLT3, 

ITRIPT(J)  , ALLOW, TURi\2,PAYTIM 

GO  TO  lOLO 

45  IF(8O0.-T0TAL)5O,46,46        .■'... 

46  PE<NLTY  =  800. -TOTAL 

47  PAYTIi-'  =  PFNLTY  +  TOTAL 

48  ALL0W=?EMLTY+PENLTl  +  PENLT2  +  TF(L)  +  Tii(L)+TF(J)+TB(  J) 

49  N'=M  +  1 

52 


5  0 
51 


1 


GC    IC,    4C 

IF (9CG.-TCTAL) 1C00»51 »51 

PFNLTY=U 

GC    TO    4  7 

SUN'  =  TRIPT(L)+TRIPTU)+TF(L)+T5(  J)+TuRN+GAP 

IF(SUN'-75C..  )  l>;0n»5n2»5]] 
50?  ALLCW  =  GAP-TB(L)-Tr  (  J)+800.-.SUN' 
5''^?  ,v  =  K  +  1 

TCTAL  =  ?;UM-TURN  '         ■  •      " 

PAYTIM  =  8Uv..  • 

5  05    ei.CCKX  =  H'LCCK(L>^<-]Or;.+BLCCK(  J) 

RC:UTFX  =  R0UT£(L)*10C.+RCUTE(  J) 

PUMCH5o9  ,  i\  .ROUTEX  »  8LGCKX  ,  TF  (  L)  .  T  lMCi\  (  L  )  >T  IMCFF  (  J  )  »  Tb  (  J  )  ,  TOTAL  >ALLC 
IWtTURNjPAYTIM 

GC  Tc  1000  "  ,;   . 

511  IFISUr-'-HOO.  )  5C2>5C2>bl2 

512  IF(Sljy-9JO.)513,513»1000 

51?  ALLCW  =  GAP-TB(L)-TF.(  J)      ,  v 
v  =  M  +  "! 

TCTAL  =  S!..iv.-TURN 
pAYTIf^=«^UM 
■  GC  TC  505         "■ 
1000  CCVTINUE  "'='  'i   •  :..^v. 

L  =  L  +  1 
IF(L-IPW) 11»11 »1005  • 

1.00  5  STOP 

8  FCR>'AT(9X,2I4,F5.1  »2r6.0»F5.1»5X,F6.0) 
41  FCRMAT(3X,2I5,F5.] »2 I5,F5. 1»7X,F6.0) 
509  FCPMAT( I3,2I5,F5.1»  2I5,F5.1,7X,   F6 .0 , 2F7. 2 , F7 .  1) 
43  FCRMAT(  I3,2I5»F5.1  >  2  I  5  »  F5  .  1  ■>  F7  .  1  .  r6.  C  .  2F7.  2  »  F7.  1) 
FND 
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PRCGRAy       HZ.     4 
TOTAL    PAY    PROGRAM 

101  OVERfO. 
PENLTY=0 

103    RFAD1»R0UT1.8L0K1 

I -^C  ROUT  1-100.  )  102  »  103  » 103 

102  Rt:AD?vV,ROUT2»RLO<2»PEMLT»PAY 
BL0X  =  R0UT1*10C.+RL0'<1 
ROUX=ROUT2*100.+BLOK2 
IF(PAY-80  0. ) 15.9,10 

10  OVrR=:(PAY-800.  )*.5 
9    Ir(PEiMLT)  15»1A,  11 

14  PFNLTY=0 
GO    TO    12 

11  PE.MLTY  =  PENLT^.5 

12  PAYX=OVER+PENLTY+PAY 
"^UMCHS.N.BLOX.ROUX.PAYX 
60    TO    101 

]     FORMAT     (3X,^I5) 

?    FORMAT(I3,?I5,?OX,I5»22X,F7.1) 

3    FORMAK I4,2I6»F7.1 ) 

15  STOP 
END 


K»»^waf.ja-^<.»  i».t_MiM 
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PRCGRAM   NO.  5 


A  SFLFCTION   PROGRAM 


DIMENSION    RLOKK  1000)  .ROUT]  (1000)  »  PAY  (1000)  »3L0,<X(75)  ,ROUTX(75) 
1YX(75) 
RFAD100,fvlL      ....,_., 
ICn    F  ^RMAT(  15  )  ■'  •'  1 

DJ22I=1 .NL   ■  ^     ■      - 

READl.ciLOKK  I  )  .ROUTl  (  I  )  ,PAY(  I  ) 

1  FORMAT (316) 
22  CONTINUE 

44.4  RFA014,N 

14  FORMAT ( 13) 
D0333  L=1»N 

MX=0  .  •         . 

CHFCK  =  U.O        [■: 
D04  J  =  l  ,NL 
^F(BL0K1(L)-BL0K1(J) )2,4,2 

2  IF{R0UT1(L)"R0UT1 ( J) )44»4,44 
44  IF(CHECK)  15»f)5.6 

6  DCb    M=1.J< 

IF(R0UT1( J)-ROUTX(M) )ll,U,Tf 
77  IF(R0UT1 ( J)-BLOKX(M)  )5,4,5 
5  CONTINUE 

D07  .M  =  1,J< 

IF{BL0K1(J)-BL0KX(M) )88,4,88 
8  8  IF(RL0'<]  {  J)-ROUTX(V)  )7,4,7 

7  CONTINUE     . 
Cits    VX=MX+1 

66    3L0KX(MX) =BLOK] (J) 

R0UTX(MX)=R0UT1( J) 

PAYX(MX )=PAY( J) 

JK  =  MX 

CHEC!<=1.0  ,.■„      , 

4    CONTINUE  ^'' 

PUNCH12.BL0i<l(L  )  » ROUTl  (L)  .PAY(L) 

IF(MX)15.15»34 
34    PUNCH12»  (BLOKX(MX)  »ROUTX(MX)  ,PAYX(|V|X)  ,   ■MX  =  1»J<) 
12    F0RMAT(?0X,3I5.45X) 

PUNCH    16  ,■-..>■,;. 

15  FORMAT!///)  •■    .    ' 
33  3  CONTINUE  •   ' 

GO  TO  444        '  ...    ..  '-' 
15  STOP 

END  •  ■    -r    • 
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PROGRAM   ,NC.  6  ■ 

A  SELECTION  OF  UNUSED  PIECES  PROGRAM 

DIMENSION  BLOK( 100) ,ROUT( 100) 

RFAD  3»N 
3  FORMAT (  13  )  '     '    -  . 

DO  5  1=1 »N  ■   . 

PFADl »RLOK(  I  )  ,ROUT(  I)   '  '  '  . 
1  FORMAT. (?nx,2 15) 

5  CONTINUE 

■  7  READ8»   ROUTE, BLOCK, TF,TIMON,TIMOFF,TR,TRIPT 

8  FORMAT ( 9X, 2  I  A, F 5.1 ,2F5.0,F5.1,5X,F6.0) 
BL0K.1  =  ROUTE*  100. +5  LOCK 

C  ^  6  J  =  1,N 
Ir(dLOi<l-BLOK(J)  )6»7,6 

6  CONTINUE 

20  D09  •<=1,N- 
IF(BL0.<1-R0UT(K)  )9»7,9 

9  CONTINUE 

21  PUNCH8,  ROUTE, BLOCK. TF,TIM0N,TIM0FF,TB,TRIPT 

GO  TO  7 
END 


.i(Jf*i-'' 
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program  mc.  7 

a  cckputaticn  of  values  program  for  eg.  2 

Dimension  spread; 20 )   '  •  ■ 
3  format { 13 )         ■  . 
33  formak f^.o)      .   ■ 

2  FORMAT {F3.0) 

4  FORMAT{3F20.1C) 

'^EAD3»'N 

J01J=1»N 

PFA03  3,SPRFmO(J) 
1  CONTINUE      v-.    .  • 
22  RFAD?»PLAT 

DO  4A  J=1,N  • 

Y =4207.75 1 12-445. 394 19*L0GF( SPREAD (J) ) +. 06529*PLAT 

A=EXPF(Y) 

PUNCH4,Y,SPREAD( J) .PLAT 
44  CONTINUE  '.  '   •  ■- 

GO  TO  22  5   • 

END 
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PROGRAM  NO.  8. 

A  CCMPUTATICN  CF  VALUES  PROGRAM  FOR  EO.  1 


DIN'ENSICN  SPRFAD(20) 
FORMAT  (  13)        .:  '    ■   _ 
FORMAT { F4.0)      ,   .  ■  •  / 
FORMAT ( F3.C> 
F0RMAT(3F20.10) 
READ3»N 

c  ■;iJ=i»N 

Ri:AD33.SPREAD(  J)  ; 

CONTINUE  ■  !  "  •   . 

READ?»PLAT 

DO    44    J=1,N 

Y=9.583  59+.C5147*L0GF(PLAT) 

A=EXPF( Y) 

PUMCH4,A,SPREAD( J) jPLAT 
CONTINUE 
GO  TO  22 


-.4122  9*L0GF(SPREAD( J) ) 


END 
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The  objective  of  this  thesis  is  to  develop  a  mathematical  model 
for  predicting  the  total  pay  to  the  operators  of  the  public  transit 
industry  in  general  functional  forms  and  parameter  values.  Economic 
measurements  were  applied  to  each  variable  and  the  degree  of  relation- 
ship was  determined.  There  are  usually  several  constraints  and 
restrictions  in  labor  management  agreements  of  public  transit  companies. 
Platform  time  and  spread  time  are  the  two  main  constraints  on  which 
controversy  is  always  going  on  between  management  and  trade  unions. 
It  is  felt  that  platform  time  and  spread  time  are  the  two  main  deter- 
ministic variables.  These  variables  are  used  to  establish  a  likelihood 
that  they  have  a  relationship  to  the  total  pay  time. 

By  conducting  a  two-way  analysis  of  variance  test  it  is  concluded 
that  on  the  total  pay  time  the  affect  of  spread  time  is  more  signifi- 
cant than  the  affect  of  platform  time.  The  computation  was  carried  on 
IBM  1620  computer  and  the  total  pay  time  functions  were  plotted  for 
each  combination  of  parameter  values. 


